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Abstract 

[J Cooperative beamforming in relay networks is considered, in which a source transmits to its 
destination with the help of a set of cooperating nodes. The source first transmits locally. The cooperating 



^ ■ nodes that receive the source signal retransmit a weighted version of it in an amplify-and-forward (AF) 

' fashion. Assuming knowledge of the second-order statistics of the channel state information, beamforming 

weights are determined so that the signal-to-noise ratio (SNR) at the destination is maximized subject 
to two different power constraints, i.e., a total (source and relay) power constraint, and individual relay 
power constraints. For the former constraint, the original problem is transformed into a problem of one 
lO ' variable, which can be solved via Newton's method. For the latter constraint, the original problem is 

in 

o. 

^D . this case, it is shown that when the number of relays does not exceed three the global solution can always 



transformed into a homogeneous quadratically constrained quadratic programming (QCQP) problem. In 



be constructed via semidefinite programming (SDP) relaxation and the matrix rank-one decomposition 
technique. For the cases in which the SDP relaxation does not generate a rank one solution, two methods 



H I are proposed to solve the problem: the first one is based on the coordinate descent method, and the second 

one transforms the QCQP problem into an infinity norm maximization problem in which a smooth finite 
norm approximation can lead to the solution using the augmented Lagrangian method. 

Index Terms 

Cooperative beamforming, channel uncertainty, relay networks, fractional programming, semidefinite 
programming. 
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I. Introduction 

Cooperative beamforming (CB), also called distributed beamforming has attracted considerable research 
interest recently, due to its potential for improving communication reliability. One form of distributed 
beamforming, the so-called distributed transmit beamforming, is a form of cooperative communications 
in which a network of multiple transmitters cooperate to transmit a common message coherently to a 
Base Station (BS). The distributed transmit beamforming can provide energy efficiency and reasonable 
directional gain for ad hoc sensor networks HI, Q. The challenges and recent progress of distributed 
transmit beamforming are discussed in ||3l. Another form of distributed beamforming is the distributed 
relay beamforming, in which a set of cooperating nodes act as a virtual antenna array and adjust their 
transmission weights to form a beam to the destination. This can result in diversity gains similar to those 
of multiple-antenna systems Q, ifTOl . Various effective cooperation schemes have been proposed in the 
literature, such as amplify-and-forward (AF), decode-and-forward (DF) lH, coded-cooperation ISJ, and 
compress-and-forward [6|. The AF protocol, due to its simplicity, is of particular interest fTOl . 

In distributed relay beamforming, the objective is to determine source power and beamforming weights 
according to some optimality criterion. Existing results for this problem can be classified into those that 
rely on channel state information (CSI) availability at the relays Q, lH, 111, and those that allow for 
channel uncertainly, i.e., that rely on statistics of CSI, such as the covariance of channel coefficients, 
or imperfect CSI feedback ifTOl . ifTTl . |[T2l . as opposed to explicit CSI. The latter class of techniques 
is particularly important because CSI is never perfectly known at the transmitter. This work picks up 
on some important results presented in 1.10.1 . in which a source transmits a signal to a destination with 
the assistance of a set of AF relay nodes In |[T0'|, the problem of obtaining the beamforming weights so 
that the signal-to-noise ratio (SNR) at the destination is maximized subject to certain power constraints 
is considered, i.e., individual relay power constraints and a total power relay constraint. For the case of 
individual relay power constraints, a semidefinite programming (SDP) relaxation plus bisection search 
technique was proposed in ifTOl . When the SDP relaxation generates a rank-one solution, then this is 
the exact solution of the original problem; otherwise, the exact solution cannot be guaranteed, and the 
authors of |[TOl proposed a Gaussian random procedure (GRP) to search for an approximate solution 
based on the SDP relaxation solution. However, GRP is time-consuming and sometimes ineffective. 

In this paper, we investigate the same scenario as in |fTOl|, i.e., cooperative beamforming under the 
assumption that the second-order statistics of the channel state information (CSI) are available. The 
beamforming weights are determined so that the SNR at the destination is maximized subject to two 
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different power constraints: (i) a total (source plus relay) power constraint, and (ii) individual relay 
constraints. The differences of this work as compared to ifTOl . are the following. 

• Our first kind of power constraint includes the source power as well as the power of the relays. In 
a wireless network all nodes have power constraints, therefore, placing a constraint on the source is 
more realistic. However, this results in a more difficult optimization problem. A similar constraint 
was also used in lfT3l . For this case, we transform the original problem into a problem of one 
variable, which can then be solved via Newton's method. 

• The second kind of power constraint is exactly the same as that of fTOl, but our work contributes 
new results and more efficient algorithms to reach the solution. In particular, 

- We show that when the number of relays does not exceed three, the global solution can always 
be constructed via SDP relaxation and the matrix rank-one decomposition technique. 

- For the case in which the SDP relaxation solution has rank greater than one, we propose two 
methods to obtain an approximate solution that is more effective than the Gaussian random 
procedure employed in lITOl . The first method is based on the coordinate descent method. The 
second method transforms the original problem into an infinity norm maximization problem, for 
which a smooth finite norm approximation results in a solution using the augmented Lagrangian 
method. 

• For both types of constraints, we obtain exact solutions for the special cases in which the channel 
coefficients between different node pairs are uncorrelated and follow a Rayleigh fading model. These 
cases were not discussed in ifTOl . 

The remainder of the paper is organized as follows. The mathematical model is introduced in ^ 
In ^im the SNR maximization subject to a total power constraint is presented. The SNR maximization 
subject to individual relay power constraints is developed in Section ^IVI Numerical results are presented 
in ^ to illustrate the proposed algorithms. Finally, ^Vll provides concluding remarks. 

A. Notation 

Upper case and lower case bold symbols denote matrices and vectors, respectively. Superscripts *, T 
and t denote respectively conjugate, transposition and conjugate transposition. | • | denotes the amplitude of 
a complex number. det(A) and Tr(A) denote determinant and trace of matrix A, respectively. Amin(A) 
and Ainax(A) denote the smallest and largest eigenvalues of A, respectively. A ^ and A ;^ mean 
that matrix A is Hermitian positive semidefinite, and positive definite, respectively. A ^ B denotes 
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that A — B is a positive semidefinite matrix. rank(A) denotes the rank of matrix A. diag(v) denotes 
a diagonal matrix with diagonal entries consisting of the elements of v. ||a|| denotes Euclidean norm 
of vector a. I„ denotes the identity matrix of order n (the subscript is dropped when the dimension is 
obvious). E(-) denotes expectation. 



Relay 1 



Desination 




Relay N 



Fig. 1. System model. 



II. System Model and Problem Statement 

The system model is the same as in ifTOJ and is depicted in Fig. [T] It consists of a source node, a 
destination node and N relay nodes, each node equipped with a single antenna. The source transmits 
signals to the destination with the help of relay nodes. We assume that the direct link between the source 
and destination is very weak and thus ignored. The channel gains from the source to the ith relay, and 
from the ith relay to the destination, are denoted respectively by fi and gi. 

Communication between source and destination occurs in two stages (slots). During the first stage, 
the source broadcasts its signal to the relays. During the second stage, the relays working in AF fashion 
transmit a weighted version of the signal that they received during the first stage. Let ^/P^ s be the source 
signal, where Pg is the source transmit power and s is the information symbol with E(|sp) = 1. The 
received signal at the ith relay is given by 

Xi = \/^s fiS + Vi (1) 

where Vi represents the noise at the ith relay having zero mean and variance a^. The ith relay weights 
the received signal and transmits Zi = WiXt where Wi is the weight. The received signal at the destination 
equals 

N N N 

y = ^giZi + iy = \/p's ^ WifigiS + ^ WigtVi + u (2) 



j=i 



j=i 



j=i 
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where v is the noise at the destination having zero mean and variance a^ . 

Let us assume that the second-order statistics of the channel gains /j's and gi% are known. We also 
assume that /j and g^, Vi,j are statistically independent. Define 

h= [/i5i,--- jNQNf, 

g = [gi,--- ,9Nf, 

R = E{hht}, 
Q = E{ggt}, 



andD = diag(E{|/i|2},... ,E{|/^|2}). 



(3) 



In general, Q and R are full matrices. In case of uncorrelated Rayleigh fading, in holds that E(/*/j) = 0, 
and K{g^gj) = 0, \/i ^ j, in which case R and Q both are diagonal. 
From ([2]), the signal component power is given by 

2"! 



P, = E 



and the total noise power Pn equals 



N 



^s^Wifigis 



i=l 



P„ = E <^ ^ WigiVi + 
The SNR at the destination is given by 



Af 



j=l 



P,wtRw 



a + a w^Qw. 






R virtRw 



Pn O-^l + Vi^tQw' 

The total relay transmit power and transmit power at the ith relay are respectively given by 



Af 



Pr = J]]E{|zip} = P^w^Dw + o-VW 



j=l 



E{\zi\^} = {P,Dii + a^)\wi 



(4) 

(5) 

(6) 

(7) 
(8) 



where Da is the (i, i)th entry of D. 

Our goal in this paper is to determine the beamforming weights Wi's such that Td is maximized subject 
to certain power constraints. In this paper, we consider two kinds of power constraints. The first kind 
corresponds to the case in which the total power of the source and all relays is constrained, i.e.. 



Ps + Pr<Po- 



(9) 
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where Pq is the maximum allowable total transmit power of the source and all relays. The second kind 
is the individual relay power constraints in which each relay node is restricted in its transmit power, i.e, 

Pr,^ < Pi (10) 

where Pi is the maximum allowable transmit power of the ith relay. 

III. SNR Maximization Under Total Power Constraint 

From ^ and Q, the SNR maximization problem subject to a total power constraint is expressed as 

Ps wtRw 

max — ^ r— — (11) 

Pa,w CT^ 1 + wTQw 

S.t. Ps + PgW^'Dw + Cr^wW < Pq. 

We give the following lemma, the proof of which can be found in Appendix lAl 
Lemma 1: Let P° be the solution of the following 

s.t. < P, < Po 
where 

Si = R-V2DR-1/2 ^ (^2/p^)j^-l^ (^3) 

and S2 = R"^/2QR"^/2 ^ {a"^ / Pq)^^^ , (14) 

Let w° be the eigenvector associated with the smallest eigenvalue of P°Si + {Po — P°)S2- Then (P°,'w°) 
is the solution to the problem of diil) . 

Remarks: Here we assume that R ^ 0. If R ^ 0, the methodology is similar. In fact, from Appendix |Al 
the problem of ([TT]) is also equivalent to 

max Ps{Po - P.)A^ax ([P^D + a^I + (Pq - P,)Q]-5R[P,D + a^I + (Pq - P.)Q]-^) (15) 

s.t. 0<Ps< Pq. 

A similar procedure can be used to solve the above problem. 

Let us normahze Pg by letting x = Pg/Po, < x < 1. With this, the problem of (fT2ll is equivalent to 

Po x(l-x) 
max — . . r^-^ (16) 

^ cr^ Amin(2;Si + (1 - X)S2) 

S.t. 0<x<l. 
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A. Si and S2 are both diagonal 

In case of uncorrelated Rayleigh fading, R and Q are diagonal matrices. Then, Si and S2 are both 
diagonal, and as it will be shown next the exact solution can be obtained analytically. 

By denoting the {k, k)-t\\ entry of Si and S2 as a^ and bk, respectively, the problem of (IT6l ) becomes 

min Amin ( :; Si + -S2 

0<a;<l y 1 — X X 

( CLk bk 
= mm mm < 1 

0<x<l k=l,-,N [ 1 — X X 

( ak bk 
= mm mm < 1 

k=l,-,N 0<a;<l [1 — X X 

= _min (^/ofe + v/ftfc)^ 

= (^/^+\/^)'- (17) 

The above minimum is attained for 

V ^0 /io\ 

X = j= (18) 



where 



B. Si or S2 is not diagonal 



^ + V ^^=0 



fco = arg min (Vofc + V^fc)- (19) 

fc=l,--- ,A'' 



Lemma 2: The optimal x of fl76l) lies in [xi,Xu] where 

XI = -^^, (20) 

1 + Vc 

Vd 

and Xu = 1= (21) 

1 + Vd 

where c = Amin(S7 828^ ' ), and d = Amax(Si' ' ^2^^ )■ 

The proof is given in Appendix |B] 

From Lemma |2l to solve the problem of (fT6l ) is equivalent to solving the problem of 

min Amin ( :; Si + -S2 ) (22) 

X \\ — X X J 

S.t. Xl < X < Xu- 

The objective in (l22l ) is in general not a convex function over [2;/,x„]. We will use Newton's method 
to search for the stationary points. Let us start by denoting 

G(x) = Si + -S2, X G [xi,Xu]. (23) 

1 — X X 
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Note that G(x) depends smoothly on x E [x^, x^] as any order derivative of G{x) exists. We assume that 
G{x) has a simple spectrum for x £ [xi, x„]. This is a reasonable assumption for general Si and S2 (see 
|[T4l . ifTSl §4]). Under this assumption, Amin(G(x)) also depends smoothly on x G [xi,Xu] IH. First- 
and second- order necessary conditions for x to be a local minimizer are respectively fTD. Theorem 2.2, 

2.3] 

— A,nin(G(x)) = 0, (24) 

and ^Ainin(G(x)) > 0. (25) 

If (I25] ) holds with strict inequality, then x is a strict local minimizer ||3T1 Theorem 2.4]. In Newton's 
method, the {k + l)th iteration is given by [311 Ch. 3] 

-^■^min[G[X)) 

Xk+1 =Xk- "fc-^f- , fe = 0, 1, • • • (26) 

gp'AminiG^Xjj 

where a^ > is chosen such that x,t+i does not exceed [x^jX^], and otherwise, a^ ■<— aA:/2- 

In the iteration expression (l26l ). we need to calculate the first- and second- order derivatives of 
'^min(G(x)). Let uo(x) be the eigenvector associated with Amin(G(x)). Let Ufc(x), /c = 1, • • • , A^ — 1 be 
the eigenvectors associated with the other eigenvalues Afc(x) of G(x), respectively, where Ai(x) > • • • > 
AAr_i(x) > Amin(G(x)). The first- and second- order derivatives of Amin(G(x)) (the so-called Hadamard 
first variation formula and Hadamard second variation formula ifTSl §4]) are respectively given by I1I61I . 

El 

AA^i,(G(x)) = uo(x)t^^uo(x), (27) 

dx dx 

H^'x fr( ^^ , ,t^!^M (, ^^ 2|u,(x)t^uo(x)f 

and ^A^..(G(x)) = uo(x)t^-^uo(x) - ^^ ;,^.(,) _ ;,^^^(g(,)) (28) 

where 

dG(x) ^ 1 g^ _ J_s^^ ^29) 

dx (1 — x)^ x^ 

d^G(x) 2 ^ 2 

and , „ = 7- t^tSi + ^82. (30) 

dx^ (1 — xj'^ x'^ 

IV. SNR Maximization Under Individual Relay Power Constraints 
From ^ and (ITOl l. the SNR maximization problem subject to individual relay power constraints is 

expressed as 

Ps wtRw 

w (T^ 1 + wTLJw 

s.t. {PsDkk + (T^)\wk\^ < Pk, k G I 
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where / = {1, 2, • • • , N}. The problem of (|3T]) belongs to the class of quadratically constrained fractional 
programs. In ifTOt . this problem was analyzed and an SDP relaxation plus bisection search technique was 
proposed. Here, we first consider the case of uncorrelated Rayleigh fading, and show that an exact solution 
can be obtained. Then, for the general fading case, we propose two methods that are more efficient than 
the search method of |10|. As it will be shown in the simulations section, the random search approach, 
in addition to being time consuming, can result in a noticeable performance gap as compared to the 
proposed approaches. 



A. R and Q are both diagonal 

By using the Dinkelbach-type method ifTSl . we introduce the following function: 

F(t) = max 

w |_ ■ a^ 

s.t. {PsDkk + (T^Mwkl'^ < Pk, kel. 
The relation between F{t) and the problem of (OTI ) is given in the following property ifTSll . 

Property 1: 

(i) F{t) is strictly decreasing, and F{t) = has a unique root, say t*; 



f{t, w) = ^w+Rw - t(l + w+Qw) 



(32) 



(ii) Let w* be the solution of i\32\l corresponding to t*. Then av* is also the solution of ( l3iD with the 
largest objective value t* exactly. 

According to Property [TJ we aim to find t* and the associated w*, which is also the solution of (|3T| |. 
To this end, by denoting the (A;, A;)th entry of R, Q as r^, qk, respectively, we rewrite 

^ fPs \ 

f(^t,w) = -t + '^i-^rk-tqkj\wk\'^ (33) 

k=i ^"^ ^ 



to get that 



associated with the optimal 



) = -t+ 


^ ft {P. , 




Ph Ps^. +„. -^ n 


Wk\^ = < 


P,Duu+<T^ a^^k tqk>U 

otherwise 



(34) 



where 



(35) 



. , X X > 
^{x) ^ { (36) 

otherwise. 
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To find the root of F{t) = 0, let us denote 

tk = ^, k = l,--- ,N (37) 

and their rearrangement 1 1 < 1 2 < • • • < i at corresponding to r^, qk, Pk, and Dkk, respectively. With 
these, we rewrite (l34l) as 

^ P / P \ 

Fit) = -t + Y: -^rr^^n ^'^ - ** ) • ^3^) 

k=l 



PsDkk + a^ V^" 



Note that F(0) > and F(iAf) = -In < 0. Thus, it follows from Property [U that < t'' < In- The root 
t* is determined based on the following theorem, the proof of which is given in Appendix [Cl 

Theorem 1: If Fitkg) = for an integer ko, then t* = tkg- Otherwise, let ko be the smallest integer 
such that F{tk„) < 0. Then 

t*= i+y ^'"^ \ y^fhE^ . (39) 

Once t* is obtained, we can obtain w* from (|35]) . 

B. R or Q is not diagonal 

1) Equivalent QCQP and SDP relaxation: The problem of (|3TI ) is equivalent (up to scaling) to a 
QCQP, as stated in the following lemma. The proof of the lemma is given in Appendix iDl 

Lemma 3: Let w° be the solution of the following homogeneous QCQP problem: 

max w^Rw (40) 

w 

s.t. w'I'AfcW < 1, A; G I 

where 

PsDkk + fr'^-, ^ ,.,. 

Afc = Jfc + Q (41) 

and Jfc is a matrix with all zero entries except for the (k, k)th entry one. Let 

PsP>kk+o- otT ° M^^ 

77 = max w JfcW . (42) 

kal Pk 



Then 4=w° is the solution to the problem of (|J7]). 

Remarks: In fact. Lemma [3] states that the QCQP of (l40l ) and the problem of (OTT ) are equivalent up to 

scaling. 
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Note that the constraint in (l40l ) is convex but the objective is concave. Thus, the problem of (l40l ) is 
not a convex problem. In fact, this problem belongs to the class of problems involving maximization of 
convex functions over a convex set [19|. 

The SDP relaxation is a popular method for QCQP problems. Let X = ww^^, and we can write 
viftRw = Tr(RX), w^A^vif = Tr(AfcX). With this, we can rewrite the problem of ( |40l ) as 

min -Ti-(RX) (43) 

s.t. Tr(AfcX) < 1, ke I 

X^O, 

rank(X) = 1. 
Dropping the non-convex constraint rank(X) = 1, we obtain the SDP relaxation ll30l 

min -Tr(RX) (44) 

.X. 

s.t. Tr(AfcX) < 1, k€ I 

X^O. 

The SDP of (l44l ) a convex problem which can be effectively solved by CVX software ll32l . Let X* be 
such a solution. Obviously, if X* has rank one, then it is the solution to the problem of (l43l) and hence 
generates the solution to the problem of (l40l ). Otherwise, a search technique may be used to obtain the 
suboptimal solution of the original problem, e.g., the Gaussian random procedure (GRP) ifTOl . For general 
R and Q, the solution X* from CVX software does not necessarily have rank one (in fact, for general 
R and Q matrices, the SDP of (l44l) does not necessarily have a rank one solution). Some examples on 
the above claim will be given in the simulation section below. 

The SDP relaxation problem of (l44l ) has several advantages as compared to the SDP relaxation of 
Uni. First, it obtains the same objective value while avoiding the bisection search. Second, for A^ = 2, 3, 
it attains the global optimal solution in polynomial time. In other words, for A^ = 2, 3, one can ensure 
that the problem of (l44l ) has a rank one solution. Moreover, one can construct a rank one solution from 
any non rank one X* in polynomial time. In fact, for A^ = 2,3, the problem has been solved using the 
complex matrix rank-one decomposition |[20l Theorem 2.1], as stated in the following theorem. 



Theorem 2: For N = 2,3, the problem of H4\i has a rank one solution. Let X* be any one of the 
solutions. If^* has a rank greater than one, one can construct a rank one solution from X* in polynomial 
time by using the complex matrix rank one decomposition. 
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For the case in which the solution X* from the CVX software has a rank greater than one, the GRP 
can used, although it is in general time-consuming and sometimes ineffective. In the following, we give 
two more effective methods for that case. 

2) Coordinate descent method: If the solution X* from CVX software has rank greater than one we 
can use the coordinate descent method [21, §8.9], |[22l §2.7], |[23l . |[24l to directly deal with the original 
problem of (|3TI ). Note that the constraints of the problem of (|3T] ) are some bounds for the elements of 
w, i.e., a Cartesian product of some closed convex sets (see |22, §2.7]). The idea behind the coordinate 
descent method is the following. At each iteration, the objective is minimized with respect to one element 
of w^ while keeping the other elements fixed. The method is particularly attractive when the subproblem 
is easy to solve (e.g., there is a closed form solution) and also satisfies certain condition for convergence 
|[22l Proposition 2.7.1], 1231 Theorem 4.1], 1241 §6]. The coordinate descent algorithm applied to our 
problem is as follows. 

Algorithm 1: 

1) Set e = 10^^; Choose an initial point wr^; Set k = 0. 

2) For p = 1 : N, determine the optimal pth element while keeping the other elements fixed. This 
results in Wp/ 

3) w'^'+i = w^; 

4) // | | ^.|| < e, stop; 

5) k = k + l; Go to 2). 

In the following, we show that the subproblem stated in Step 2 has a closed form solution (see Theorem 
[3l) and also study its convergence to a stationary point (see Theorem |4l). 

It is easy to verify that minimizing the objective with respect to the kth element of w while keeping 
the other elements fixed leads to the following optimization problem: 

ailvl"^ + hy + bly* + ci 

y 02|yp + 622/ + &2y* + C2 

s.t. |y| < /3 



max , , ^ , ^^ ^ , (45) 



where /3 = yJPk/[PsDkk + o"^), ai = Rkk, «2 = Qkk and hi,h2,ci,C2 can be inferred from (OTI ). For 
example, when k = l,\et v^ = [y, w^]^ and 



Q= ^^^ ' , andR= ' . (46) 
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Then 61 = wHi, ci = w"fQiw, 62 = ^^^2 and C2 = 1 + w"fRiw. 

For the solution of (|45] | we give the following theorem, the proof of which can be found in Appendix 

E 

Theorem 3: If ai/a2 = 61/62 = ci/c2, the objective in ( 1451 ) /^ a constant, and the optimum y, i.e., 
y*, is any value satisfying \y\ < /3. Otherwise: If the equation (oi — ta2)/3^ + 2|6i — ^62!/? + ci — tc2 = 
has a real root, i.e., ti, such that \hi — tib2\ > (^1^2 — oi)A then the optimal y is given by 

/ = (3e-''' (47) 

where 61 G (— 7r,7r] is the argument of bi—tib2; Else, let ^2 be the rootof\bi—th2\^ = {cLi—ta2){ci—tc2) 
such that 1 61 — ^2^2 1 < (^2^2 — ai)/5, then the optimal y is given by 

^.._\bi-t2b2\^.;o^ (48) 



^202 — ai 

where 62 is the argument of 61 — ^262- 
Remarks: The roots ti and t2 in Theorem |3] can both be obtained in closed form. 

For the coordinate descent method, obviously the function value sequence converges. However, in 
general additional conditions for convergence to a stationary point (or fixed point used in |24, §6]) are 
needed. 

Theorem 4: The sequence {w } generated by Algorithm Ul converges globally to a stationary point. 

Proof: Our proof is based on (2T. Proposition 2.7.1] and its proof. Let us denote the objective in 
(I3T]) as /{w). Let w = {wi, ■ ■ ■ , wn) be the limit point of the sequence {w*^}. We first show 

/(w) > f{wi,W2, • • • , wn), Vwi. (49) 



If f{wi,W2, • • • , Wn) is a constant, then obviously (1491 ) holds. If f{wi,W2, ■ ■ ■ , wn) is not a constant, 

to see why, let us assume that ( |49l ) does not hold. A verbatim repetition of the proof for Il22l Proposition 
2.7.1] results in 

/(w) = f{wi + ewi,u)2,--- ,wn), Ve G [0,eo] (50) 

for some vi / 0, eo > 0. But from Theorem [3l (l50l ) does not hold for any t;! / 0, eo > if 
f{wi,W2, ■ ■ ■ ,wn) is not a constant. Thus, ( |49l ) holds. Similarly, we show 

/(w) > /(tDi,--- ,iBj-i,Wj,iI)j+i,--- ,wn), ywj (51) 

for J = 1, • • • , A'^. This completes the proof. ■ 
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3) p-norm approximation: If the solution X* from the CVX software has rank greater than one, we 
can also use j»-norm approximation plus an augmented Lagrangian method to solve the problem of (l40l ). 
The convergence of the augmented Lagrangian method can be found in |[3T| . First, we can show that the 
problem of (l40l) is equivalent (up to scaling) to 



mm 

w 



max w^Ai.w 



(52) 



s.t. w^Rw = 1. 



To see why this is the case, let w* be the solution to the problem of (1521 ) associated with the optimal 
objective value max^g/ w^^^A^w* = C. Then, ^w* is the solution to the problem of (l40l ) associated 
with optimal objective value (-^w*)'l"R(-^w*) = ^. Otherwise, let us assume that the solution to 
the problem of (l40l) is w' with w'^Rw' = C'l > i. Thus, w" = ^==w' satisfies vi^'^Rw' = 1 and 



maxfcg/ w"^AkMv" = -^ < C. This contradicts the optimality of w* for the problem of (l52l ). In fact, 
the two problems are equivalent up to scaling. 



On denoting Di = diag( J^^^^±^, • • • ^J^^^fr^), u = Diw, Ri = D^^^RD^^ and Qi = 
D^i^^QDj^^, we rewrite the problem of ( [52l ) as 

min u^Qiu + ||u||^ (53) 

u 

s.t. u' Riu = 1 

where ||u||oo = max^g/ |njt| is the infinity norm. Note that ||u||oo is not smooth |[26ll . However, we can 
approximate ||u||oo by (smooth) p-norm, i.e., ||u||p = (X]fce/ Kfcl^)^''^' ^o that ETll . ESl 

||u||oo = lim ||u|L, (54) 

and ||u||oo < ||u||p < iVVP||u||oo. (55) 



When p is sufficiently large, the approximation is good. In fact, from (I55l ). it is easy to show that given 
a tolerance e, the relative error does not exceed e as long as p > log A^/log(l + e). For example, for 
AT = 10, e = 1%, we get p > 232; for N = A<d, e = 0.5%, we get p > 740. 

Now, using ||u||| , p > 1 as a smooth approximation to ||u||^, we turn to solve the following 

min u^Qiu + ||u||L (56) 

s.t. u^Riu = 1. 



We use the augmented Lagrangian method BTl §17] to solve the problem of ( |56l ). Since the augmented 
Lagrangian method was originally proposed for real variables, we first modify our problem as follows. 
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Define ESI 

' Re(u) \ 

, (57) 

Im(u) / 

Re(Qi) -Im(Qi^ . 
Im(Qi) Re(Qi) 

, Re(Ri) -lm(Ri) , 
K = I I , (59) 

Im(Ri) Re(Ri^ 

Jfc \ 
and Jfc = I (60) 

Jfc / 
where Jfc is defined in Lemma [3l and Re(-), Im(-) denote the real and imaginary part respectively, then 

u^Qiu = z^Fz, (61) 

u^Riu = z^Kz, (62) 

and ||u||2p = ( ^(z^Jfcz)M . (63) 

With these, we rewrite the problem of ( |56l ) as 

min z'^Fz + (/)p(z) (64) 



s.t. z^Kz -1 = 



where (/>p(z) is defined as the right hand side of (l63l) . 

Now we can apply the augmented Lagrangian method, given by 

L(z; A; ii) = z^Fz + (/.p(z) - A(z^Kz - 1) + — (z^Kz - 1)^ (65) 



where A is the Lagrangian multiplier, and the fourth term in the right hand side of (1651) is the penalty 
function. The algorithm is described as follows: 

1) Choose an initial estimate X^^' of A* and fi = 0.001. Set k = 1. 

2) Determine Zfc to be a minimizer of L(z; A'^^^^; h); 

3) Compute A('=) = \^^-^^ - {^lY.Zk - 1)//^; 
A) If a convergence test is satisfied, stop; 

5) k = k + l; Go to 2). 
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For Step 2, we use the backtracking line search Newton's method with Hessian modification IIBTI 
Algorithm 3.2]. The iteration expression is 

^(i+i) = ^{i) + api (66) 

with Pi = -(V^L + /3I)-^VL (67) 

where /3 is chosen such that V^L + /3I is positive definite, e.g., /3 = Amin(V^L) + 10~^, and a is the 
step size determined by the backtracking line search described as follows 131] Algorithm 3.1]: 

a) Set a = 1, ci = 10"^ p = 0.5; 

b) Repeat: if L{z^^^ + apf, •; •) > L(z(*); •; •) + ciapf VL, then a ^ pa. 

In the algorithm, we need to calculate VL and V^L given by 

2 
VL = 2Fz + V0„-2AKz + -(z^Kz-l)Kz (68) 

and V^L = 2F + v2(/>p - 2AK + -(z^Kz - 1)K + -Kzz^K. (69) 

The calculation of Vcl)p and V'^ipp is given in Appendix El 

Remarks: For the initial estimate A^*^^ of A*, note that when p = 1, then (f)p{z) = z^z and A* can be 

expressed in closed form as Amin(K^^/^FK~^/^ + K^^). We choose this as A*^"-*. 

V. Numerical Results 

In this section, we provide some examples illustrating the proposed algorithms. For more simulation 
results on beamforming itself the reader can refer to lITOl . We consider a channel model as follows: 



h = h + ViHfi (70) 

and gj = gj + ^/^ gj (71) 

where fi and gj are means, tpi and ipj are variances, fi and gj both are zero-mean random variables 
with unit variance. We assume that fi, fj, gi and gj, \/i ^ j are independent, /i = corresponds to the 
scenario in which there is no line-of-sight (LOS) path (Rayleigh fading), while fy^O corresponds the 
scenario in which there is an LOS path (Rician fading). Thus, the matrices D, R and Q are given by 

D = diag(|/i|2 + Vi,---,l/7v|' + V'^), 

Qij = 9i9*j + y/WPj^ij', 



and Rij = {fif* + ^J^iMj 5ij){gig* + ^^pupjSij) 



where 5ij is the Kronecker function. 
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A. SNR maximization under total power constraint 

Please refer to ^III] for details. First, we consider a network consisting of A^ = 6 relays with channel 
parameters given by 

/i = 0.2202 + 0.8130i, /2 = -0.4075 - 0.7644i, 
h = -2.0107 + 0.40161, U = -0.4503 + 0.0678i, 
/s = 0.8588 - O.llSOi, h = -0.1219 + 0.4260i; 

ipi = 3.8042, V'2 = 2.6326, ^3 = 4.7590, 
tpi = 0.4989, V's = 1.2576, ^g = 1-2484; 

g^ = -0.3726 + 0.8007i, 52 = 0.4592 - 0.20451, 
^3 = -0.8769 + 0.46711, 54 = -0.9270 + 0.54301, 
55 = -0.0063 - 0.49771, % = -0.7783 - 0.77121; 

ipi = 0.3913, ip2 = 0.4791, y^s = 0.0865, 
(^4 = 2.7813, (^5 = 4.8960, 996 = 4.6789. 

Fig. |2]plots Amm(G(x)) for x in [2;i,Xu] = [0.1711,0.7077] with 100 uniform points. Using Newton's 
method for starting points xq = xi, xq = Xu, the convergent points (0.2156,1.2191), (0.5844,1.2694) 
are also plotted in Fig. |2] Fig. |3] plots the iteration process under the stopping test: F'°+i~^'° I < 10^'^ 
and |^Ainin(G(x))| < 10~^. It can be seen from Fig. [3] that Newton's method converges rapidly. 
Second, we consider a network consisting of A^ = 6 relays with channel parameters given by 

/i = -0.4751 + 0.73401, /2 = -0.0449 - 0.46091, 
h = 0.0239 - 1.51541, /4 = 0.5130 - 0.17551, 
/s = -0.2017 + 0.67171, h = 1-0134 - 0.19851; 

V'l = 2.4707, V2 = 3.9193, V's = 2.4121, 
^A = 3.8879, ^5 = 1.2050, -06 = 3.0901; 

51 = 0.5360 - 1.29321, 52 = 1-7471 - 0.89141, 
53 = 0.0955 - 0.15771, 54 = -0.6795 + 0.24791, 
55 = 0.5815 + 0.50391, qq = -0.3090 + 0.84131; 

ifi = 3.9655, ^2 = 0.2693, tp^ = 0.9205, 
(^4 = 0.5567, c^5 = 3.3901, ipQ = 2.9367. 

Fig- HI plots Amin(G(x)) for x in [xi,Xu] = [0.2754,0.6392] with 100 uniform points. Using Newton's 
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method for starting points xq = xi, xq = Xu, the same convergent point (0.4087, 0.6060) is also plotted in 
Fig.H Fig.|5]plots the iteration process under the stopping test: |^^^*^^^| < 10"^ and \^Xmm{G{x))\ < 
10^^. It can be seen from Fig. [5] that Newton's method converges rapidly. 



B. SNR maximization under individual relay power constraints 

Please refer to ^TV] for details. In ifTOl . the authors stated that, based on their simulations, the SDP 
relaxation always has a rank one solution. However, no analytic proof was provided for that claim. 
However, although a rank one solution often occurs, for general R and Q the SDP relaxation does not 
necessarily have a rank one solution. This can be seen in the following examples, for which the SDP 
relaxation has a rank greater than one. 

First, we consider a network consisting of A^ = 4 relays with 

/ 



Q = 



and R = 



V 



2.1 


.73 + .751 


.43 + 1.11 


.70 - .331 ^ 


73 - .751 


1.6 - 


-.20 + .181 


.57 -.711 


43-1.11 - 


-.20 -.181 


2 - 


-.52 -.451 


70 + .331 


.57 +.711 - 


-.52 + .451 


.98 j 


1.6 


-.74- .161 


.084- .571 


-.19 +.671 ^ 


-.74 +.161 


1.1 


-.88 +.311 


-.44- .241 


.084 + .571 


-.88 -.311 


2 


.20- .141 


-.19 -.671 


-.44 +.241 


.20 +.141 


1.5 j 



(72) 



(73) 



For simplicity, we let Di = I (defined in ^IV-B3I ) and denote the SDP relaxation solution from CVX 
software by X*. The eigenvalues of X* are 

0.0000, 0.0000, 0.2064, 1.8148. 

Thus, X* has rank two rather than rank one and can be eigen-decomposed as 0.2064uiuj + 1.8148u2U2 
where ui and U2 are eigenvectors associated with the eigenvalues 0.2064 and 1.8148 respectively. We 
obtain the objective values of the problem of (l40l ) for SDP relaxation, GRP from ifTOl . coordinate descent 
method from ^IV-B2[ and p-norm approximation from ^IV-B3I (starting points: \/l-8148u2 or some 
samples from (JAA(0, X*)) as, respectively: 

SDP relaxation: 3.74112 
GRP (10*5 samples from CAA(0,X*)|- 3.6970 
Coordinate descent method: 3.7076 
p-norm approximation: 3.7069 (p = 1024J 
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It can be seen that the objective values from GRP, coordinate descent method and p-norm approximation 
are close to each other (with a difference < 0.3%) and close to the SDP relaxation solution (with a 
difference < 2%). It can be seen that: although the GRP attains a close performance compared with 
the other two methods, it is time consuming in the sense that it needs much more time (processes 10^ 
samples from CJ\f{0, X*)). The augmented Lagrangian L(z; A; fi) (defined in (l65l )) during the iteration 
is plotted in Fig. [6l The objective value during the iteration for the coordinate descent method is plotted 
in Fig. |7] It can be seen that for these two algorithms the iteration converges rapidly. 
Second, we consider a network consisting of A^ = 6 relays with 
/ 



Q 



.778 -.658 -.6461 .135 +.2691 -.273 + .0051 

-.658 + .6461 2.20 -.379-1.141 .253 - .8721 

.135 -.2691 -.379+1.141 2. .689 + .2981 -.547 

-.273 -.0051 .253 +.8721 .689 - .2981 1. -.655 

.088 + .2611 -.337-1.021 -.547 +.1601 -.655 -.1921 

-.021 + .0131 .444 +.0351 .373 - .6931 .132 + .1071 -.721 



and 



R 



3.44 -.263 +.0541 .572+1.731 .490 - .2761 

-.263 -.0541 3.09 -.342-1.491 .926+1.131 

.572-1.731 -.342+1.491 2.70 -.493 + .8651 

.490 +.2761 .926-1.131 -.493 -.8651 3.09 

-.613 + 1.621 -.282 +.7131 -.396 -.8261 .541 - .3301 

-.014 -.3751 -.211 -.9111 .149 +.8361 -.552 + .2211 



088- .2611 -.021- .0131 \ 

337+1.021 .444 -.0351 

373 + .6931 

132- .1071 

721 - .2761 



.1601 
.1921 
2.40 
.2761 



1.09 



(74) 



.613- 1.621 -.014 +.3751 \ 

.282 -.7131 -.211 + .9111 
.396 +.8261 .149 -.8361 

.541 + .3301 -.552 -.2211 

2.75 -.442 - .3521 
.442 +.3521 2.08 



(75) 



The eigenvalues of X* are 

0.0000, 0.0000, 0.0000, 0.0000, 0.8369, 2.3774. 

Thus, X* has rank two rather than rank one and can be eigen-decomposed as 0.8369uiu| + 2.3774u2U2 
where ui and U2 are eigenvectors associated with the eigenvalues 0.8369 and 2.3774 respectively. We 
obtain the objective values of the problem of (l40l ) for SDP relaxation, GRP, coordinate descent method, and 



p-norm approximation (starting points: a/2.3774 U2 or some samples from CAf{0,X.*)) as, respectively: 

SDP relaxation: 9.33816 

GRP (l{f samples from CJ\f{^,X*)): 8.1472 
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Coordinate descent method: 8.9428 
p-norm approximation: 8.9409 (p = 1024J 

It can be seen that: the objective value from GRP has a significant (> 10%) difference from the 
SDP relaxation solution; p-norm approximation and coordinate descent method attain objective values 
close to each other; the improvement of objective value from coordinate descent method is (8.9428 — 
8.1472)/8.1472 = 9.77% compared with GRP. It can be seen that for this example, GRP is time 
consuming and ineffective in the sense that it needs more time (processes 10® samples from CJ\f{0, X*)) 
but attains worse performance compared with the other two algorithms. The augmented Lagrangian 
L(z; A; /i) (defined in ( [65l) ) during the iteration is plotted in Fig. [8] The objective value during the 
iteration for the coordinate descent method is plotted in Fig. |9] We can see that for these two algorithms 
the iteration converges rapidly. 

VI. Conclusion 

We have investigated the problem of cooperative beamforming under the assumption that the second- 
order statistics of the channel state information (CSI) are available. Beamforming weights are determined 
so that the SNR at the destination is maximized subject to two kinds of power constraints. The first 
kind of power constraint is a constraint on the total power, i.e., source plus relay power. The second 
kind of power constraint is a constraint on each relay's transmit power. For uncorrelated Rayleigh fading 
scenario, we attained the exact solution. For generic fading scenario, we focused on the case in which 
the SDP relaxation does not produce a rank-one solution and proposed two methods to solve it. The 
numerical simulations suggest that the proposed methods are more effective than the method of ifTOl . 

Appendix A 
Proof of Lemma [T] 

Let (P°, w°) be the solution to the problem of (HB. We can show that P° + P°w°''Uw° + a'^w°''w° = 
Pq. Otherwise, let us assume that P° +P°w°'f'Dw° +a'^w°'<w° < Pq. Let /3 = {PQ-P°)/{P°-iv°^T>^N° + 
cr^w°W°), and hence /3 > 1. It is easy to verify that (P°, y^w") satisfies the constraint but results 
in a larger objective value. This violates the optimality of (P°,w°). With this, the problem of ([TT]) is 
equivalent to 

Ps wtRw 

max — ^ r— — (76) 

Pa,w cr^ 1 + wTQw^ 

s.t. Ps + Psw"^Dw + cr^virtw = Po- 
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It follows from the constraint in (1761) that 



wt(P,D + o-2l)w 



By using (TTTI ). we rewrite the problem of (1761) as 

Ps (Pq - P,)wtRw 

max —77 — TT ?; ; r^ri — ('o) 

P.,w a2wt[P,D + cj2l+(Po-i^s)Q]w 

s.t. Ps + PgW^'Dw + o-^wW = Pq. 

Note that the objective in ( TTSl ) has the same value at w and /3iw, V/3i 7^ 0, w 7^ 0. Thus, the problem 

of (TTSI ) is equivalent to 

P, (Pq - P,)wtRw 

™^w a2wt[P,D + cT2l + (Po-P,)Q]w ^ ^ 



s.t. < Ps < Po, w / 0. 



,2n 



Ct2. _ „,a2 



(T^I = Ps— I + (Po-P.)— I, (80) 

Po Po 



Further, we rewrite 



which enables us to write 

P,D + a^I + (Po - P,)Q = P, Td + ^ 1 + (Po - P,) I Q + ^ 1 . (81) 

With this, by using the fact that for Ci ^ and C2 ^ 1291 p. 549] 

1 xtCix 

max — j— - — , (82) 



a^\ .„ „./„ ^2 



A,,i„(C-^/'C2C-'/') -^0 xtC2x' 
the problem of ( 1791) is equivalent to the problem of (lT2l) . 



Appendix B 
Proof of Lemma [2l 

Obviously, neither x = nor x = 1 is the solution to the problem of ( [T6l ). Define the function 

xSi + (l-x)S2 Si S2 ^/n1^ /oa^ 

-^W = n ^ = -^ + —, x£ (0, 1 . (83) 

X[l — X) 1 — X X 

Let X G (0, 1) and Ax 7^ be an increment such that x + Ax € (0, 1). Using Taylor series expansion, 

we approximate 

1 1 Ax (Ax)2 

+ T. ^ + .. .\. , (84) 



l-(x + Ax) 1-x (l-x)2 (1-6)^' 
1 _ 1 Ax (Ax)2 
X + Ax X x2 ^: 



and ^ = — H -3 — (85) 
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where ^i and ^2 both he between x and x + Ax. From (1831 ). (1841 ) and (1851 ). we get 

K{x + Ax) = ir(x) + Ax ( (^-^Si - ^82) + (Ax)2 ( (T^^i + ^82). (86) 



Note that the third term in the right hand side of (I86b is positive definite. By using the facts that Il29[ p. 
549] 

cSi < S2, and dSi >z S2 (87) 



it is not difficult to prove that 



1 1 ^0 xG (0,Xi] 

Si - ^82 <^ - (88) 



{'^-xf x2 I ^0 xG [x„,l). 

With these, we know that: if x G (0, x;] and Ax < 0, then K{x + Ax) >- K{x) and it follows from 
Weyl's inequality |33, p. 181] that Amin(-^(a^+Ax)) > \^in{K{x)); if x G [x„, 1) and Ax > 0, similarly, 
^ram{K{x + Ax)) > \^iji{K{x)). This completes the proof. 

Appendix C 
Proof of Theorem [H 

If ko = 1, then < t* < ti, and 

%fc-i''gfc>0, k = l,--- ,N. (89) 



Thus, F(t) = in (|38]) leads to 



^ A 



k=l ^ s kk 

The desired result can be obtained from the above equation. 
If ko > 1, then 4o-i < t* < tk^, and 



t* + y ^ -^fk - t*qk = 0. (90) 



Ps . > k = ko,--- ,N 

^h-t*qk { (91) 



^" ' < A: = I,--- ,A;o-l 

Thus, F(i) = in (I38]l leads to 

^ Pk fP. 



k=ko 

The desired result can be obtained from the above equation. 
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Appendix D 
Proof of Lemma [3] 

Note that the constraints in (|3T]) can be rewritten as 



Pk 



-W 



^JfcW < 1, k £ I. 



(93) 



Note that the objective in ( [3T] ) has a greater value at aw than that at w, Va > 1, w 7^ 0. Thus, there 
exists j £ I such that ({PsDjj + a"^) / Pj)w'^ J jW = 1, i.e., at least one constraint is active. With this, 
the constraints in (|3T]) can be rewritten as 

PsDkk + T^ 



max 
fee/ 



s^kk 

Pk 



-W^JlW = 1. 



By using (|94l ). we rewrite the problem of ( [3T] ) as 

P, wtRw 



max 



w (T^ maxfcg/ wt AfeW 
s.t. max ^ — ^^^^— ■w^JfcW = 1. 



(94) 



(95) 



PsDkk + cr^ t ■ 
fee/ Pfc 



Note that the objective in 
of (I95] ) is equivalent to 



has the same value at w and /3iw, V/3i / 0, w / 0. Thus, the problem 
P, wtRw 



max 



w a^ maxfcg/ wt A^w 
s.t. w/0. 



(96) 



Similarly, the problem of ( |96l ) is equivalent to 



max — ^w'Rw 



w a^ 



(97) 



s.t. maxw^Ai,w = 1. 
fee/ 



Obviously, the problem of ( 1971 ) is equivalent to the problem of (140] ). This completes the proof. 



Appendix E 



Calculation of V^p and V'^4>p 



We have 



V0P 



9: 



^^E 



fee/ 



0p(z) 



p-i 



JfeZ 



(98) 
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and 

.2. _ d 



V 



^ \^ </.p(z) J 0p(z) \dz J \dz J 0p(z) ^ \^ 0p(z) y 

(99) 

Appendix F 
Proof of Theorem [3] 

By using the Dinkelbach-type method lITSl (cf. ^IV-AK we introduce the function 

F(t)=max /(t,y) (100) 

y 

s.t. \y\ </3 

where 

fit, y) = ai|y|2 + biy + h\y* + ci - t{a2\y\^ + 62^ + bly* + cs). (101) 

Similarly to Property [T] in ^IV-A[ F{t) is a strictly decreasing function and the equation F{t) = has a 
unique root t*. The optimal y* associated with F{t*) is also the solution for the problem of (1451 ) with 
the optimal objective value t*. 

To obtain the expression of F{t), we denote y = |y|e'^ and write 



id I tl, 4-1, \* ^ — 'iS 



fit, y) = (ai - ta2)\yV + (&i - th2)e'' + (61 - t62)*e 



|y| +ci -tc2 



< (ai - ta2)\yr + 2|6i - t62||y| + ci - tc2. (102) 

The equality in (11021 ) occurs when the argument of bi — tb2 equals —9 (if 61 — ^62 = 0, then 9 is arbitrary). 
With this, we let r = |y| and write 

F{t) =ci - tc2 + max (ai - ta2)r'^ + 2\hi - tb2\r (103) 

r 

S.t. < r < /3. 

Further, it is easy to get: 

1) When ai — ta2 > 0, i.e., t < ai/a2, the optimal r is /3 (unique), and we wirte 

F{t) = (ai - ta2)P^ + 2|6i - tb2\^ + ci - ^2. (104) 
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2) When ai — ta2 < 0, i.e., t > ai/a2, the optimal r is given by (unique) 

r* = minji^l^M, ^1 (105) 

[ to2 - ai J 

and we write: If t > ai/a2 and |&i — ^62! > (i«2 — ai)/?> then 

F(t) = (ai - ta2)/3^ + 2\bi - tb2\f3 + ci - tc2; (106) 

If i > ai/a2 and |6i — t62| < (^02 — ai)/3, then 

|6i -t52p 

^W = ^^ ^ + ci-tc2. (107) 

1:02 — ai 

3) When oi — ta2 = 0, i.e., t = 01/02, we know: If 61 — ^62 / (i.e., 61/62 / 01/02), the optimal r is 
/3 (unique), and 

F(ai/a2) = 2|6i - (ai/a2)62|/3 + ci - (01/03)02; (108) 

If 61 — ^62 = (i.e., 61/62 = 01/02), the optimal r is arbitrary in [0, /3], and 

i^(ai/02) = Ci-(0i/02)c2. (109) 

Recall that F{t) is a strictly decreasing function and the equation F{t) = has a unique root i*. Thus, 
one and only one of the equations (fT04l ). (IT06l ). (ITOTl) . (IT08] ). ( fT09l ) satisfies F(t*) = 0. 
Based on the analysis above, it is not difficult to obtain the desired result. 
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Fig. 2. Amin(G(a-)) for x in [xi, Xu] = [0.2754, 0.6392] with 100 uniform points; The left point is for starting point xo = xi 
and the right point is for the starting point xo = Xu', SNR = 10 dB; A total power constraint. 




Fig. 3. The iteration process; The upper line is for starting point xo = xi and the lower line is for the starting point xo = Xu', 
SNR = 10 dB; A total power constraint. 
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Fig. 4. Amin(G(a;)) for x in [a::;,a::„] = [0.2754, 0.6392] with 100 uniform points; The same convergent point is for starting 
point xo = xi and xo — Xu\ SNR = 10 dB; A total power constraint. 
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Fig. 5. The iteration process; The upper and lower lines converge to the same point for starting point a;o = xi and a;o = x^, 
SNR = 10 dB; A total power constraint. 
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Fig. 6. The augmented Lagrangian L(z; A; /i) (defined in J65t ) during the iteration process of the proposed algorithm; Individual 
relay power constraints; p-norm approximation. 
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Fig. 7. The objective value during the iteration process of the proposed algorithm; Individual relay power constraints; Coordinate 
descent method. 
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Fig. 8. The augmented Lagrangian L(z; A; jj.) (defined in J65j) during the iteration process of the proposed algorithm; Individual 
relay power constraints; p-norm approximation. 
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Fig. 9. The objective value during the iteration process of the proposed algorithm; Individual relay power constraints; Coordinate 
descent method. 



June 1, 2010 



DRAFT 



